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We analyze the effect of spin degree of freedom on fidelity decay and entropy production of 
a many-particle fermionic(bosonic) system in a mean-field, quenched by a random two-body in¬ 
teraction preserving many-particle spin S. The system Hamiltonian is represented by embedded 
Gaussian orthogonal ensemble (EGOE) of random matrices (for time-reversal and rotationally in¬ 
variant systems) with one plus two-body interactions preserving S for fermions/bosons. EGOE are 
paradigmatic models to study the dynamical transition from integrability to chaos in interacting 
many-body quantum systems. A simple general picture, in which the variances of the eigenvalue 
density play a central role, is obtained for describing the short-time dynamics of fidelity decay and 
entropy production. Using some approximations, an EGOE formula for the time {tsat) for the onset 
of saturation of entropy, is also derived. These analytical EGOE results are in good agreement with 
numerical calculations. Moreover, both fermion and boson systems show significant spin dependence 
on the relaxation dynamics of the fidelity and entropy. 
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I. INTRODUCTION 

Investigation of non-equilibrium dynamics and sta¬ 
tistical relaxation in isolated, interacting, many-body 
quantum systems has emerged as an important research 
area in the last few decades [1-3]. Magnificent exper¬ 
imental progress, particularly in the field of ultracold 
atoms, has made it possible to artificially realize inter¬ 
acting many-body quantum systems to understand var¬ 
ious quantum phenomenon. Systematic understanding 
of non-equilibrium dynamics of such complex systems is 
still not available. This is not only a fundamental require¬ 
ment but may also be crucial for applications in quan¬ 
tum information [Ij. A complex many-particle quan¬ 
tum system exhibits non-equilibrium dynamics following 
a change in one of the system parameters, usually referred 
to as quantum quench. Coherent quench dynamics has 
been observed for both bosonic [4] and fermionic [5] sys¬ 
tems. The main questions addressed in these studies are 
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to understand: (i) whether and how a quantum system 
thermalizes, (ii) universal aspects of the dynamics, (iii) 
the relation between localization, integrability and ther¬ 
malization, and (iv) behavior of various observables after 
equilibration. Although theoretical advancements have 
been made to understand these questions using spin lat¬ 
tice models [6-16], deriving universal properties for such 
complex systems remains an open question. 

It is now generally accepted that statistical relaxation 
in the quench dynamics of an isolated interacting many- 
body quantum system is closely related to quantum chaos 
[17-19]. Signatures of quantum chaos appear in the sta¬ 
tistical behavior of eigenvalues and eigenfunctions [20]. 
Quantum integrable systems have been found to relax 
to an equilibrium state characterized by the generalized 
Gibbs ensemble [10, 12-16, 21, 22]. On the other hand, 
non-integrable quantum systems relax towards thermal 
equilibrium. Therefore, for these systems, for exam¬ 
ple, expectation values of various few-body observables 
approach their long time average values given by the 
Gibbs ensemble. The eigenstate thermalization hypothe¬ 
sis (ETH) is considered to be the underlying mechanism 
for thermalization in isolated quantum systems [17, 18]. 
ETH states that the eigenstate expectation values of typ¬ 
ical observables are smooth functions of energy eigenval¬ 
ues with off-diagonal elements being exponentially small 
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in system size. Various aspects related to ETH have been 
studied by many groups using lattice models of interact¬ 
ing spins, fermionic as well as bosonic; see for example 
[11, 23-27], 

In order to understand the role played by quantum 
chaos in attaining thermalization of isolated finite quan¬ 
tum systems, thermodynamical description of isolated 
interacting quantum systems using random interactions 
was first addressed in [28]. EGOE are paradigmatic mod¬ 
els to study many-body quantum chaos [29]. The con¬ 
stituents of real physical systems interact predominantly 
via two-body interactions in the presence of an aver¬ 
age one-body (mean-field) interaction. Flambaum and 
Izrailev started investigation of time evolution of generic 
quantum many-body systems [30, 31]. Going further, 
EGOE(l-|-2) generated by one plus two-body interactions 
was employed in [32] to investigate thermalization in iso¬ 
lated quantum systems using ergodicity principle. One 
of the most significant aspects of EGOE(l-|-2) is that the 
system undergoes a transition from regular to chaotic re¬ 
gion with increasing strength of the random two-body in¬ 
teraction quench that deeply affects eigenvalue and eigen¬ 
function structure [33]. These random matrix ensembles 
are used so far to study (i) ergodicity principle for expec¬ 
tation values of several different operators(observables); 
(ii) determine a region of thermalization using the cri¬ 
terion of equivalence between different definitions of en¬ 
tropy and temperature; (iii) representability of occupan¬ 
cies of single particle states using Fermi-Dirac (or Bose- 
Einstein) distribution; (iv) calculation of expectation val¬ 
ues using canonical distribution. See [32, 34, 35] and 
references therein. It was also shown that the onset of 
Wigner’s surmise in spectral fluctuations is not sufficient 
for finite systems to thermalize [32]. 

EGOE are generic, though analytically difficult to deal 
with, compared to lattice spin models as the latter are 
associated with spatial coordinates (nearest and next- 
nearest neighbor interactions) only. Moreover, the uni¬ 
versal properties derived using EGOE extend easily to 
systems represented by lattice spin models, for exam¬ 
ple eigenvalue density is Gaussian in the mean-field basis 
unlike a classical GOE [11]. Due to two-body character 
of the interactions, the three different types of non-zero 
many-particle matrix elements (zero, one and two par¬ 
ticle transfer) can be derived analytically by action of 
creation and annihilation operators on the many-particle 
states in the Fock-space representation with the system 
Hamiltonian (see Eq. (1) ahead) expressed in second 
quantized form. These non-zero matrix elements will be 
a linear combination of two-body matrix elements (cho¬ 
sen to be GOE) and hence, correlated. Due to this, ma¬ 
trix structure for an EGOE is different from a classical 
GOE which accounts for many-body interactions among 
system constituents (although fluctuation properties of 
EGOE are same as that of a GOE) [33, 36]. 

In addition to particle number, spin (S') quantum num¬ 
ber is important for realistic systems like atomic nu¬ 
clei, quantum dots, nano-metallic grains, ultracold spinor 


gases etc. As group symmetries define various quan¬ 
tum numbers, in general, one has to consider EGOE 
with group symmetries. The first non-trivial and im¬ 
portant (from the point of view of its applications) 
embedded ensembles are EGOE(l-|-2)-S [BEGOE(l-|-2)- 
E] with spin degree of freedom for fermions [bosons] 
[37, 38]. In the notation BEGOE(l-|-2)-E, ‘B’ stands 
for bosons and F is their fictitious spin-i degree of free¬ 
dom. EGOE(l-|-2)-S' has been successfully used to an¬ 
alyze universal conductance fluctuations in mesoscopic 
systems [39]. Also, BEGOE(l-|-2)-F may be useful in 
exploring general structures of spinor condensates [40]. 

We analyze unitary evolution of an isolated many-body 
quantum system followed by an interaction quench. Ini¬ 
tially, the system is prepared in one of the eigenstate 
of unperturbed mean-field Hamiltonian h{l). Dynam¬ 
ics starts with instantaneous (in a time interval much 
shorter than any characteristic time scale of the model) 
change of interaction strength A from zero. Then, the fi¬ 
nal perturbed Hamiltonian is H = h{l) + XV{2) with 
eigenvalues E; V(2) is the random perturbation. To 
characterize system evolution, we analyze fidelity decay 
and information entropy for small times. There have 
been several experimental [41] and analytical [42] studies 
showing that fidelity/survival probability exhibits non¬ 
exponential behavior at long times. Fidelity is related 
to Loschmidt echo and is defined as the overlap between 
the initial state and its evolved counterpart. Shannon 
entropy is a measure of complexity that captures partic¬ 
ipation of other (^ initial) states in system dynamics. 
Fidelity decay and entropy production after a quantum 
quench have been considered in a number of works; see 
for example[23, 31, 43-48]. Going further, we analyti¬ 
cally and numerically analyze the effect of many-particle 
spin on fidelity decay and entropy production. We also 
derive analytical formulas for saturation value and sat¬ 
uration time for entropy. Estimation of saturation time 
is important for development of algorithms for quantum 
optimal control [49]. Now we will give a preview. 

Section 2 briefly introduces the EGOE(l-|-2)-S' and 
BEGOE(l-|-2)-E models along with the chaos markers 
generated by these ensembles. Spin dependence of spec¬ 
tral variances is also investigated in detail as it explains 
the spin dependence of fidelity and entropy. Section 3 
gives definitions and analytical formulas for fidelity, in¬ 
formation entropy, saturation value and saturation time 
for information entropy. Numerical results are compared 
with analytics in Section 4. Finally, Section 5 gives con¬ 
clusions. 


II. EMBEDDED FERMIONIC/BOSONIC 
ENSEMBLES WITH SPIN 

To analyze time dynamics of an isolated system 
initially prepared in a mean-field basis state (unper¬ 
turbed) followed by a random two-body interaction (non- 
integrable) quench, we model the system Hamiltonian 
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by EG0E(l+2)-S' for fermions [and BEGOE(l+2)-E for 
bosons]. Following our previous works [37, 38] , we briefly 
introduce EGOE(l+2)-S' [also BEGOE(l+2)-F] gener¬ 
ated by one plus two-body random interactions preserv¬ 
ing total spin S. 


A. Definition and construction 

We consider a system of m (m > 2) fermions dis¬ 
tributed in n number of single particle (sp) orbitals each 
with spin s = ^ so that number of sp states is N = 212. 
Distributing these m fermions in N sp states gener¬ 
ates a complete set of basis states spanning the Hilbert 
space. As s = i, the two-particle spin s = 0 or 1 and 
many-particle spin S = m/2, m/2 — 1, ..., 0 or 1/2. 
EGOE(l-|-2)-S' is defined by the Hamiltonian H that pre¬ 
serves total m particle spin S, 

H = h{l) + XV{2) ; H(2) = W=°(2)0H"=i(2) . (1) 

Here, A is the strength of the two-body interaction V{2). 
We choose one-body Hamiltonian h{l) = nj, spec¬ 
ified by sp energies Ci with unit average level spacing. 
Note that, rii are number operators acting on the sp 
states f = 1, 2, ..., 12. As two-particle spins can 
take two values, the two-body Hamiltonian V{2) is a 
direct sum matrix of matrices in spin 0 and 1 spaces 
with dimensions 12(12 -|- l)/2 and 12(12 — l)/2 respec¬ 
tively. Thus, V(2) is defined by two-body matrix ele¬ 
ments that are independent of the mg quantum number. 
Then the m particle states can be classified according to 
U{N) D 17(12) G SU{2) algebra with SU{2) generating 
spin S. As H preserves S, it is a scalar in spin SU{2) 
space. 

EGOE(l-|-2)-S' is generated by the action of H on 
many-particle basis states, with H®=°(2) and V‘^^^{2) 
being independent GOEs in two-particle spaces. Given 
the sp energies and the two-body matrix elements, the 
many-particle Hamiltonian H is first constructed in the 
Ms = M™*" [M™*" = 0(i) for m even(odd)] representa¬ 
tion using the spin-less EGOE formalism and then states 
with a given S are projected using a spin projection oper¬ 
ator [37, 38]. Thus, the many-particle Hamiltonian H 
will be a block diagonal matrix with each diagonal block 
corresponding to a given spin S [37]. Similarly, one can 
define BEGOE(l-|-2)-F with s replaced by f, s replaced 
by / and S replaced by F [38]. 

From now on, whenever we discuss the results that 
apply to both EGOE(l-|-2)-S' and BEGOE(l-|-2)-F, we 
refer to them as EE(l-|-2)-5' results. 


B. Transition or chaos markers 

As a function of the interaction strength A, the em¬ 
bedded ensembles with one plus two-body interactions 
exhibit three transition or chaos markers (Ac,A^,A(i): 


(a) as the two-body interaction is turned on, level fluc¬ 
tuations exhibit a transition from Poisson to GOE at 
A = Ac; (b) with further increase in A, the strength func¬ 
tions (also known as the local density of states) make a 
transition from Breit-Wigner (BW) to Gaussian form at 
A = A^ > Ac; and (c) beyond A = A^, there is a region 
around A = Ad where basis dependent thermodynamic 
quantities like entropy behave alike. For further details of 
the three chaos markers for EE(l-|-2) for spinless systems, 
see [32, 33, 50]. It was established that EE(l-|-2)-5' also 
exhibits these chaos markers and their spin dependence 
was also deduced [37, 38]. Exact formula for the ensem¬ 
ble averaged variances [second moment of the eigenvalue 
density with fixed-(f2, m, S')] was also derived using trace 
propagation method. Importantly, the variances explain 
the spin dependence of various chaos markers. 


C. Spin dependence of spectral variances 


Initially, the system is in eigenstate 'k(2 = 0) = jfc) of 
unperturbed mean-held Hamiltonian h{l). The dynam¬ 
ics begin with the sudden change in the parameter A, 
the strength of the perturbation, in a time interval much 
shorter than any characteristic time scale of the model. 
With the quench V (2) of strength A, the state unitarily 
changes after time 2 to (2) given by 

'I'(2) = lfc(2)) = exp—*772 \k) . (2) 

Here, 77 is given by Eq. (1). Then, the probability that 
the state \k) changes to the state ]/) is H4_>/(2), 

= \{f\km^ = \Ak^f{t)f ; 

13) 

E 

where the coefficients C///_ are the overlaps of mean-field 
basis states (Jfc) and ]/)) with the eigenstates \E) of the 
post-quench Hamiltonian 77. We use 7. = 1 so that 2 is in 
E~^ units. It is important to mention that, although not 
explicitly included in the notation, the energy eigenstates 
\E) of 77 depend on many-particle spin S{F) and thus, 
Wk^f{t) and too. 

As we will see ahead in Section 3, the ensemble av¬ 
eraged spectral variance S) of the initial basis 

states jfc) is the single most crucial parameter in the 
EGOE(l-|-2)-S' [also BEGOE(I-|-2)-F] formalism to ana¬ 
lytically understand the spin dependence of fidelity and 
entropy. Firstly as argued in [51], S) will be essen¬ 

tially independent of fc and therefore for EGOE(l-|-2)-S' 
for fermionic systems, we have 


cr^(m,5) = A 


2 ^V(2) 


{m,S) 


'H 


{m,S) 


( 4 ) 


Note that ‘overline’ denotes ensemble average and 
tTy( 2 )('m, S) and a‘/j{m, S) are variances for the perturba¬ 
tion V (2) and the Hamiltonian 77 respectively. Assuming 
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that h{l) is fixed, the ensemble averaged spectral vari¬ 
ance S) follows simply from Eq. (1), 




Now, for a uniform sp spectrum having unit average level 
spacing, the follows easily from Eq. (7) of 

[37], 


CT^(i)(m, 5') = — m(S7-|-2)(r2—m/2) —2flS'(S'-|-l) . (6) 


Similarly, the ensemble averaged variance generated by 
the perturbation V(2) is = P{VL,m,S). For 

BEGOE(l-|-2)-E, Eq. (5) with S replaced by F ap¬ 
plies and the variance propagator is F). Ex¬ 

plicit formulas, derived using trace propagation method, 
for the variance propagators P{n,m,S) and Q{n,m,F) 
are given in [37] and [38] respectively. For sake of com¬ 
pleteness, we give the exact formulas for these in the 
Appendix. Again, for a uniform sp spectrum with unit 
mean level spacing, E) follows from Eq. (6) by 

transformation and applying modulo operation, 


^h{i) 


{m, F) = - 2){n + m/2) -b 20F(F -b 1) 


(7) 

For simplicity of notation, from now on, we drop the ‘bar’ 


over cr^ as we will only deal with ensemble averages. Note 
that variances are always positive as for given (n, m), the 
many-particle spin can only take the allowed values m/2, 
m/2 — 1, ..., 0 or i. 

General behavior of cr^(m, S) with S for fermionic sys¬ 
tems, using Eqs. (4)-(6) and the formula for variance 
propagator P{Q,m, S) given by Eq. (26), is shown for 
some values of (f2,m. A) in Fig. la. It is seen that 
a^(m, S) decreases with increasing spin as the number 
of states connected with each other decrease. Similarly, 
the results for BEGOE(l-b2)-E are shown in Fig. lb. As 
seen from the figure, the trend is opposite for bosonic 
systems. 

In order to understand the scaling of cr^(m, S) with sys¬ 
tem size m, we examine the asymptotic form of cr^(m, S). 
For EGOE(I-b2)-S', for fixed density po = m/fl with 
m —>■ oo and S = ml2 (the maximum value) we have. 


Similarly for minimum spin S' = 0 (assuming m is even) 
we have. 

For S varying from 0 to m/2 for m even (i to m/2 for 
m odd), the cr'^irn, S) varies between the above two lim¬ 
its. It should be noted that for an arbitrary and fixed 


S, the expansions are complicated as one has to include 
0(I/m^) terms with p increasing with S. 

This behavior is reflected in the fidelity decay and 
hence in longer time for relaxation. These are discussed 
further in the next Section. 


III. FIDELITY AND ENTROPY: DEEINITIONS 
AND EE(l+2)-S FORMULAS 

Both fidelity and information entropy are basis de¬ 
pendent quantities. In our construction of the FI ma¬ 
trices, the chosen basis states are eigenstates of both 
h{l) and S'^ operators (although all the basis states have 
Ms = we drop everywhere). 


A. Fidelity 


Probability to find the system still in the initial state 
|fc) after time t corresponds to fidelity (also referred to 
as the return or the survival probability). It is given by, 


Wk^kit) 




2 


^ exp -iEt 


FjT’^{E) exp —iEt dE 


( 10 ) 

Here E/”’‘^(E) = |C'|'p p"*’'®(E) are the strength func¬ 
tions or the local density of states and /^'"’'^(E) is the nor¬ 
malized eigenvalue density with fixed-(H, m, S). Strength 
functions give the spreading of the basis states over the 
eigenstates and thus, fidelity is the Fourier transform in 
energy of the strength function FjT’^iE). Our focus is 
on the general scenario with single peaked energy distri¬ 
bution of the initial states. 

Following our discussion in Sec. 2.2, we know that 
EE(l-b2)-S' exhibits three chaos markers {Xc, Xi, Xd)- As 
a function of the perturbation strength A, the strength 
functions are d-functions at the /i(l) eigenvalues for very 
small A. With increasing A, the eigenstates start spread¬ 
ing over all the basis states leading to delocalization in 
the Fock basis. First, they become BW in shape and fi¬ 
nally, acquire Gaussian form around and beyond X = X( 
[37]. For A > Ac, the level and strength fluctuations fol¬ 
low GOE and hence, in this region one can replace, to a 
good approximation, Fk{E) by its smoothed (ensemble 
averaged) form. As the smoothed form of Fj^{E) changes 
from BW to Gaussian with increasing A, there are four 
situations: (i) small ‘F limit where one can apply per¬ 
turbation theory; (ii) BW region; (iii) Gaussian region; 
(iv) region intermediate to BW and Gaussian forms. Re¬ 
sults for (i)-(iii) are already discussed by Flambaum and 
Izrailev [30, 31] and briefly they are as follows. 

For very small times, we have exp—iFlt ~ 
[exp—i/i(l)t] [exp—iAU(2)t]. Using ek{m,S) = 
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FIG. 1. Variation in the ensemble averaged spectral variances of the initial states al{m,S) for various (fl, m, A) values for (a) 
EGOE(l+2)-5 and (b) BEGOE(l+2)-F’ as a function of many-particle spin S{F). See text for further details. 


{k \ H \ k) ~ (fe I h{l) I k), al{m,S) = (k \ H'^ \ k) — 
e1{m,S) ~ (fc I [AF(2)]^ | fe) and the dominant terms 
(up to second order) in the expansion in time of 
exp[—iAF(2)t], the fidelity is, Wk^k{t) = 1 —S') 
The long time behavior of fidelity decay depends also 


on the shape of the distribution of F^’^{E). For A 
not far from Ac, the Fj^’^{E) has BW distribution 
with level and strength fluctuations following GOE. 
Thus, fidelity, as seen from Eq. (10), follows ex¬ 
ponential law: Wk^k{t) —f exp—r(m, S)t. Here, 
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r(TO, S') is the spreading width of the BW distribu¬ 
tion. Similarly, in the Gaussian region, we obtain 
Ak^k{t) = exp - [iek{m, S) t + <7l{m, S) /2] using 
Eq. (10). Therefore, fidelity follows Gaussian law: 
Wk^k{t) —exp—(T^(m, S) Note that when t is 
measured in [anim, S)]~^ units, the spreading widths 
r(TO, S) and crfc(m, S), respectively of the BW and Gaus¬ 
sian distributions, will be in units of aH{ni,S). In sum¬ 
mary, 


In the BW to Gaussian transition region, as demon¬ 
strated by Angom et al [52], it is possible to approximate 
Fk{E) by Student’s t-distribution in terms of the shape 
parameter a and scale parameter /3. With the transfor¬ 
mations a = -I- l)/2 and [E — ek{rn, S')] = x, 

the F^’^{E) in the BW to Gaussian transition region 
transforms to F^’^(x : v) where, 


Wk^k{t) BW_r^gion exp-r(m, S)t, 


( 11 ) 


Therefore, the logarithm of fidelity decay (lnII4_j.fc(t)) 
shows a linear dependence on time in the BW re¬ 
gion while the dependence is quadratic in the Gaussian 
regime. 


77?m,5 
^ k 


{x : v) 



(^+1) 

^ dx . 


( 12 ) 

To analytically compute the fidelity in the BW to Gaus¬ 
sian transition region, we need the Fourier transform of 
FjT’^ {x : ly). This is a topic of interest in many statistics 
problems. Using the result given in [53], we get. 


Wk-s.k(i) 


transition region 


1 (^) Jo 


t' = 


/3(i^ + l) 


2v 


t . 


(13) 


We correctly recover the results for BW and Gaussian 
limits using Eq. (13) as well: (a) for z/ = 1, we get BW 
form for Fj^’^{E) with /3 = r^/4 and (b) for i/ —>■ oo, 
we get Gaussian form with al{m,S) = /3/2. In gen¬ 
eral, <j‘f.{m,S) is related to the parameters (a,/3) by 
al{m,S) = 2a-o ^'' ^ ^ (^^) 

in [33] to analyze fidelity decay in the framework of 
EE(l-|-2). Thus, we have full EE(l-|-2)-S' theory for fi¬ 
delity decay for different regimes defined by perturbation 
strength A in terms of S dependent spectral widths. 

B. Entropy 

Information entropy S{t) (also known as Shannon en¬ 
tropy) is a measure of complexity or disorder. Assuming 
that there are total D+1 mean-field basis states such that 
/ = 0,1,..., with ]/ = 0) = jfe), the initial state in 
which the system is prepared at time t = 0, information 
entropy after time t is 

D 

Sit) = - ^ Wk^fit)lnWk^fit) . (14) 

f=o 

It is important to recognize that 

D 

'£Wk^fit) = l. (15) 

f=o 


Using Eq. (3), we have 

W.^fit) - \Cf\" 

E 

+ 2 ^ C^Cfcfcf cosiE-E')t 

E>E' 

= wz^^it) + wl^^it). 

(16) 

As system at time t = 0 is prepared in initial state jfc), for 
simplicity of notation, henceforth we denote Wk^kit) = 
Wo{t). Note that formulas for Wo{t) are derived in Sec. 
3.1 fully taking into account both the terms in Eq. (16). 

Isolated interacting many-body quantum systems equi¬ 
librate in a probabilistic way. After a long time [t —^ oo), 
the local observables fluctuate around their infinite time- 
averages with the size of the fluctuations decreasing with 
increasing system size. Therefore, it is plausible to argue 
that the second term in Eq. (16) that gives fluctuations 
for / 7 ^ 0 approaches zero for very long times. Simi¬ 
larly, in the short time limit (t close to zero), S'(t) —>■ 

al{m,S)t'^ - i-e- entropy S{t) 

will be quadratic in t. Here, Fl^^ = | (/ | | 0) are like 

transition strengths. For t neither very small and very 
large, EE(l-|-2) theory for S{t) is not available. Earlier, 
using a cascade model as an approximation, Flambaum 
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and Izrailev [31] showed that entropy initially increases 
linearly with t and eventually saturates. In this model, 
Given a set of states, these are separated into different 
classes with states directly coupled with those from the 
previous class populating the successive class. However, 
one needs to extend the Cascade model to account for 
low connectivities as well as different connectivities for 
different classes. 

To derive approximate EE(l+2)-S' formulas for S{t), 
its saturation value and the time tsat which marks the 
onset of its saturation, assume that in the sum in Eq. 
(14), variation of Wk^f{t) with / [i.e. the fluctuations in 
Wk^f{t)] is negligible. This is physical as for sufficiently 
large interaction strength A > Ac, the system becomes 
chaotic and there is no preference to any particular mean- 
field basis state. Then, one can replace H4_>/ by its 
average value W. Suppose that there are n number of 
mean-field basis states / with f ^ 0 which contribute 
to the sum in Eq. (14). Applying Eq. (15) using the 

I 

5(t) - [l 


It is quite clear that n = with S'(oo) being the 

long time saturation value of entropy. Ideally one needs 
to derive a formula for S{oo) by taking the limit t —>■ oo 
in Eq. (14). However, this is not solved yet. it is plausi¬ 
ble to expect that n should be proportional to the num¬ 
ber of mean-field basis states over which the eigenstate 
spreads. This is generally referred to as number of prin¬ 
cipal components (NPC). NPC is expected to be large 
for chaotic systems (Gaussian regime) and hence, these 
systems will thermalize unlike regular ones. The max¬ 
imal value of NPC for EE(l-|-2)-S' in Gaussian regime 
is given by putting A = 0 in Eq. (5) of [51]. There¬ 
fore, approximating n by the maximal value of NPC for 


J 


2 [I-C2(m,,5)] 


[1-C 


identity ^I = n, we have 

S(t) = -Wo(t) \nWo(t)-y"W\nW ; W= ——^ 

n 

r—l 

= -Woit)lnWoit) - [I - ITo(t)]ln ~ . 

(17) 

Now, we will derive the EE(I-|-2)-S' formulas for tsat in 
Gaussian, BW and BW to Gaussian transition regimes 
using Eq. (17). 


1. Gaussian regime 

In the Gaussian regime A > A^, fidelity Wo{t) follows 
Gaussian law; see Eq. (11). Substituting this in Eq. (17), 
we have 


— e 




In 


1 _ 


(18) 


EE(l-k2)-5, 

n ^ K(NPCinax) = Vl - ^^(m, S) ; 

3 (19) 

C^(m,5') = 1 - (Tl{m,S) . 

Here, ()^(m. S') is a correlation coefficient and k is a pa¬ 
rameter. Eqs. (18) and (19) give a theoretical description 
of time evolution of the entropy for EE(l-|-2)-S in the 
Gaussian regime. Eor very short times, expanding the 
exponential and retaining only dominant contributions 
in Eq. (18), entropy shows a quadratic growth in time. 

We define the saturation time as the time taken by the 
information entropy to reach its saturation value. Upon 
saturation, the time variation of entropy vanishes. Thus, 
setting time derivative of S{t) given by Eq. (18) to zero, 
we obtain 


^(m, S)] -I- In 


1 _ g-[l-C"(m.S)]P] 


= 0 . 


( 20 ) 


Simplifying Eq. (20), we get 


As A —oo, C^(m, S) approaches 0 and we get a lower 
bound for tsat, 


isat — 


ln(l -I- n) 

1 - C(m,S) 


( 21 ) 


j.min 

^sat 



d(m, S) N 
3 / ■ 


( 22 ) 




















m m 


FIG. 2. Asymptotic behavior of S{oo) (left panel) and tsat (right panel) with system size m. Results are shown for EGOE(l+2)- 
S ensemble with po = m/Q. = 0.5, A = 0.3 and many-particle spins S = 0 and m/2. 


Importantly, approaches 1 for A —>■ 0 giving tgat —t oo. 
Thus, within the EE(l+2)-5' formulation given here, an 
integrable system will never thermalize. As discussed in 
Sec. 2.2, an isolated closed quantum system thermalizes 
around the transition marker A = A^. In this case, the 
spreadings produced by the mean-field basis and the in¬ 
teraction basis are equal and hence, C'^{rn, S) = 0.5 that 
gives the upper bound for the time a chaotic system takes 

to thermalize, Uh — ln(At- ^™’‘^^ v^). 


2. BW regime 


In the BW region Ac < A < fidelity decay Wo{t) is 
exponential in time; see Eq. (11). Using this in Eq. (17) 
and proceeding in the same way as in Sec. 3.2.1, we 
obtain 


^ ln(l -I- n) 

^ r(m,S') 


(23) 


3. BW to Gaussian transition region 

In the BW to Gaussian transition region, Wo(t) is given 
by Eq. (13) and the formula for NPCmax was derived 
in [52]. Again, approximating n by NPCmax, we have 


n ~ K(NPCmax) 


d{m, S) 


T^{a) 


2« - 3 r2(„ _ 1) y'CHrn, ,S)[1 - CHm, ,5)] 


-.U 


13^ (2a - 3)[1 - C^(m,S')] 

2’2 2C^im,S) 


(24) 


In Eq. (24), U(-) is the hypergeometric-17 func¬ 

tion [54] and r(— ) are the Gamma functions. Time 
derivative of entropy shows that Gat should be inversely 

proportional to^^^ (1-1-^)*^. Also, x and y should be 

chosen such that tgat converges to Eqs. (21) and (23) in 
the V ^ oo and v = \ limits respectively. Using these, 


Got in the BW to Gaussian transition region is 


ign.f. — 


-fl-t-- 

[ln(l-brt)]2 V ^ 

1 ^ 3/2 


2 V tz 


;i/ = 2a-l. (25) 
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For EGOE(l+2)-S', we have also derived tedious ana¬ 
lytical expressions for S'(oo) and tgat in the asymptotic 
limit with fixed density po = to/O with m —> oo. Eig. 2 
shows the asymptotic variation of saturation values S{oo) 
(left panel) and saturation times tsat (right panel) of in¬ 
formation entropy with system size m. With increasing 
system size, both S'(oo) and tsat increase. 

Earlier, equilibration of isolated quantum system has 
been studied in the context of spin-1/2 lattice models 
conserving Sz, total spin in the z direction [11]. A re¬ 
laxation time tfi has been defined as the time required 
by hdelity to reach its long time average value. For a 
two-body interaction and initial state energy close to the 
middle of the spectrum, a lower limit for the tn was also 
estimated in [11] and this turns out to be very close to 
tsat given by Eq. (21). However, it should be clear that 
tsat is not a bound on the saturation time and it is exact 
within the EE(l-|-2)-S' formalism. Rather it turns out be 
a relatively higher estimation of the saturation time as 
can be seen from numerical results ahead. 


IV. NUMERICAL RESULTS 

Eor the numerical calculations reported for fidelity and 
entropy in this Section, we make use of the following pa¬ 
rameters. We choose fixed sp energies Ci = i + 1/i as in 
our previous studies [37, 38] and generate 20 member en¬ 
sembles. For EGOE(l-|-2)-S' ensemble with H = to = 8 
and spins S = 0 — 2, dimensions are d{m, S) = 1764, 2352 
and 720. Note that (2S'-|-1) d{m,S) = (^). For 

this choice of parameters, (a) Ac = 0.028, 0.034 and 0.05; 
(b) Xe = 0.15, 0.16 and 0.19; and (c) Xd = 0.21, 0.22 and 
0.24; for S' = 0 - 2 [37]. Similarly, for BEGOE(l-b2)- 
F ensemble with H = 4, to = 10, dimensions are 
d{m,F) = 750, 594 and 286 for F = 2, 4 and 5 respec¬ 
tively; 1) d{m,F) = Here, the 

Ac values are 0.039, 0.0315 and 0.0275 for F = 0, 2 and 
5 respectively [38]. 

First of all, the basis state energies ek{m,S) are the 
diagonal elements of the F[ matrix in the many-particle 
Fock-space basis giving efe (to, S) = (fc l/i(l) -I- AU(2)j fc). 
Note that the centroids of the ek{m,S) energies are the 
same as that of the eigenvalue (F) spectra but their 
widths are different. For each member of the ensemble, 
energies E and ek{m, S) are zero centered and scaled by 
the spectrum width a Him, S). Then, for each member at 
a given time t, \Ak^f\‘^ are summed over the basis states 
jfc) and ]/) in the energy windows ek ± and e/ ± ^ 2 - 
Here, (5i and <52 are the bin-sizes and they are chosen 
such that there are approximately 4-7 states within each 
bin. In principle, they can be different or same depend¬ 
ing on the choice of initial and final Hilbert spaces. We 
choose Si = $2 = 0.01 and ek{m,S) = 0. Numerical cal¬ 
culations, though not reported, were also performed with 
6 i = 62 = 0.05 and 0.1 and it was confirmed that the re¬ 
sults remain unaffected. Then, ensemble averaged transi¬ 
tion probability Wk^f for a fixed initial mean-field basis 


state k is obtained by binning. Note that this procedure 
is different from the one adopted in [31] where no ensem¬ 
ble average over jfc) states has been carried out. Physi¬ 
cally, the basis states indices, as used in [31], do not carry 
any significant information. However, the basis state en¬ 
ergies ekim,S) give the location of the corresponding 
strength functions and hence are meaningful [36]. Sim¬ 
ilarly, information entropy is computed numerically us¬ 
ing Eq. (14). In order to be able to compare the fidelity 
and entropy results for different spins, time t should have 
the same unit for all spins. As aH{m,S) shows signifi¬ 
cant spin dependence, we express time in units of 
with ^Ivg = [Eidim,I)\-^Y.idim,I)(j'jj{m,I)-, I \s S 
for EGOE(l-b2)-S' and F for BEGOE(l-b2)-F. 


A. Fidelity 

Variation in fidelity Wo(t) with time is shown in Fig. 3 
as a function of many-particle spin S for an EGOE(l-|-2)- 
S ensemble. Results are shown for two choices of inter¬ 
action strength: (a) A = 0.3 and (b) A = 0.1 which corre¬ 
spond respectively to the Gaussian and BW to Gaussian 
transition region for the chosen set of parameters (H, to). 
Numerical results in the figures are compared with the 
theoretical formulas given by Eqs. (11) and (13) respec¬ 
tively in the Gaussian and BW to Gaussian transition 
regimes. Agreement between numerical and theoretical 
EE(l-|-2)-S' results is excellent for short times. Inset to 
Fig. 3(b) also shows the fidelity decay for longer times. 
Fidelity decay is slower with increasing S which is in con¬ 
firmation with variation of <j\im, S) with many-particle 
spin for fermionic systems. The best-fit value of the pa¬ 
rameter z/ in Eq. (13) is found to be = 3.4, 3 and 2.6 
for S = 0, 1, and 2 respectively and the corresponding 
variances crf.im, S) are given in Table I. 

Similarly, fidelity decay for BEGOE(l-|-2)-F is shown 
in Eig. 4. Here, A = 0.2 and 0.05 correspond to Gaussian 
and BW to Gaussian transition regimes respectively. It 
is seen that agreement between numerical and analyti¬ 
cal results is excellent. Unlike fermions, fidelity decay is 
faster with increasing many-particle spin F. This is again 
in conhrmation with variation of a^{m,F) with many- 
particle spin for bosonic systems. The best-ht value of 
the parameter u in Eq. (13) is found to be z/ = 5.4, 6.6, 
and 10.0 for F = 2, 4, and 5 respectively and the corre¬ 
sponding variances cr^im, S) are given in Table I. 

As the value of the interaction strength A increase, 
the eigenstates of FI spread over all the mean-field basis 
states leading to delocalization in Eock basis. Therefore, 
hdelity decays faster with increasing quench strength A 
both for fermionic and bosonic systems as seen from 
Figs. 3 and 4. Moreover, for a given (Hjto), bosonic 
systems display faster hdelity decay in comparison to 
fermionic systems as the dimension of the fermionic 
Hilbert space is smaller due to Pauli’s exclusion principle. 
Thus, EE(l-|-2)-S' formulation captures all the essential 
features of time dynamics of hdelity decay as a function 
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FIG. 3. Fidelity decay Wo{t) with time as a function of many-particle spins S for an EGOE(l+2)-S ensemble in the: (a) 
Gaussian region (A = 0.3) and (b) BW to Gaussian transition region (A = 0.1). The points marked by red ’+’ (5 = 0), blue 
’open square’ (5 = 1) and black ’open circle’ (5 = 2) correspond to the numerical results. The solid continuous curves represent 
the corresponding theoretical results for different S values given by Eq. (11) and Eq. (13) for panels (a) and (b) respectively. 
Inset to panel (b) gives long time variation (up to t = 8) in fidelity decay in the BW to Gaussian regime. 


TABLE I. Ensemble averaged cr^ for EE(l+2)-5. 


EGOE(l-b2)-5 

BEGOE(l-b2)-E 

A 

at{m,S) 

A 

o-k{m,F) 

O 

II 

5 = 1 

CM 

II 

CM 

II 

F = 4 

F = 5 

0.1 

0.19 

0.18 

0.15 

0.05 

0.23 

0.28 

0.33 

0.21 

0.53 

0.51 

0.46 

0.06 

0.31 

0.38 

0.39 

0.24 

0.59 

0.58 

0.53 

0.1 

0.56 

0.61 

0.61 

0.3 

0.70 

0.69 

0.64 

0.2 

0.87 

0.87 

0.89 


of many-particle spin. 

B. Entropy 

Variation in information entropy with time as a fuirc- 
tion of mairy-particle spin for an EGOE(l-|-2)-5' ensemble 
is shown in Eig. 5. We choose four values of the inter¬ 
action strength: A = 0.1 (BW to Gaussian transition 
region), 0.21, 0.24 and 0.3 (Gaussian region). In all the 
cases, entropy shows a quadratic behavior for very small 
times, then switches to an essentially liirear growth aird 
eventually, saturates as predicted by Eq. (18). Fig. 5 iirdi- 
cates a good agreemeirt between numerical and analytical 


results (derived in Secs. 3.2.1 and 3.2.3). To deterirrine 
the validity of approximatmg exp[S'(oo)] by NPG, theo¬ 
retical predictions are fitted with the numerical results 
for entropy. Good agreement is found between the theo¬ 
retical predictions using n and the numerical results with 
AC = 2 in the Gaussian region and ac = 3, 3.3 and 3.6 for 
5 = 0 — 2 respectively, in the BW to Gaussian transition 
region. 

As seen from the figure, the saturation value of entropy 
for 5 = 1 is larger compared to that for 5 = 2 for a 
given interaction strength. This can also be understood 
from Eq. (18) as follows. Note that n is a function of 
C^(m, 5) and d{m,S). For a given A, d{m,S) has large 
variation with S while C‘^{m,S) is almost same for all 
spins. Therefore, the saturation value for 5 = 1 is larger 
than for 5 = 2 irrespective of A. On the other hand, 
with a given spin 5, ('^(m, 5) increases with increasing 
A leading to enhancement in saturation value of entropy 
with increasing A. 

Using the calculated value of C^{m, 5) in Eq. (19) [Eq. 
(24)] and then applying Eq. (21) [Eq. (25)], we can es¬ 
timate the saturation time tsat for different spins in the 
Gaussian [BW to Gaussian transition] region. Table II 
gives the calculated values of tsat for different 5 values. 
It is seen that tsat increases with increasing spin and this 
is consistent with the results in Fig. 5. 
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FIG. 4. Fidelity decay Wait) with time as a function of many-particle spins F for a BEG0E(l+2)-F’ ensemble in the: (a) 
Gaussian region (A = 0.2) and (b) BW to Gaussian transition region (A = 0.05). The points marked by red {F = 2), 
blue ’open square’ {F = 4) and black ’open circle’ {F = 5) correspond to the numerical results. The solid continuous curves 
represent the corresponding theoretical results for different F values given by Eq. (11) and Eq. (13) for panels (a) and (b) 
respectively. Inset to panel (b) gives long time variation (up to t = 5) in fidelity decay in the BW to Gaussian regime. 


Similarly, Fig. 6 shows time dynamics of information 
entropy for a BEGOE(l+2)-F ensemble as a function of 
many-particle spin F. We choose four different values 
of interaction strengths: A = 0.1, 0.2 (Gaussian region) 
and A = 0.05 and 0.06 (BW to Gaussian transition re¬ 
gion). In all the cases, entropy shows a quadratic behav¬ 
ior for very small times, then switches to an essentially 
linear growth and eventually, saturates as predicted by 
Eq. (18). Numerical results show good agreement with 
analytical results derived in Secs. 3.2.1 and 3.2.3 respec¬ 
tively in the Gaussian and BW to Gaussian transition 
regimes. Good agreement is found between the theoreti¬ 
cal formulas and the numerical results with k = 2 in the 
Gaussian region and k = 2.8, 2.55 and 2.3 for E = 2, 
4 and 5 respectively, in the BW to Gaussian transition 
region. Deviations from the theoretical curves are max¬ 
imum for F = 5 for all values of A and this may be due 
to the fact that dimension of the corresponding Hilbert 
space is smallest compared to other spin values. Values 
of tsat are given in Table II. The saturation times increase 
with increasing many-particle spin, just as for fermionic 
systems. 

As the parameter n is expected to depend on the nature 
of spreading of the eigenstates, it should have different 
values in the Gaussian and BW to Gaussian transition 
region. This is also in confirmation with our numerical 


TABLE II. Saturation time tsat (in (Javg units) for EE(l-|-2)- 
S'. 


EGOE(l-k2)-S 

BEGOE(H-2)-E 

A 

isat 

A 

isat 

o 

II 

CO 

S = 1 

CM 

II 

Co 

F = 2 

E = 4 

E = 5 

0.1 

7.97 

9.62 

11.49 

0.05 

6.74 

4.62 

3.28 

0.21 

3.32 

3.73 

4.21 

0.06 

5.57 

3.87 

3.15 

0.24 

3.13 

3.5 

3.97 

0.1 

3.66 

2.68 

2.2 

0.3 

2.7 

3.2 

3.65 

0.2 

3.07 

2.1 

1.7 


results. In the Gaussian regime, Eqs. (18) and (19) with 
K = 2 explain the time variation in entropy with many- 
particle spin. However, k shows dependence on many- 
particle spin in the BW to Gaussian transition regime. A 
weak dependence on spin is expected as S) varies 

with spin. The exact dependence of parameter k on A 
and many-particle spin is not explored any further in 
this paper. Thus, EE(l-|-2)-5' formulation captures all 
the essential features of time dynamics of entropy as a 
function of many-particle spin. 
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FIG. 5. Entropy S{t) as function of time t for different spins for EGOE(l+2)-S. The points, marked by red ’+’ (5 = 0), 
blue ’open square’ {S — 1) and black ’open circle’ {S — 2) correspond to numerical results. Similarly, the continuous curves 
represent the corresponding theoretical predictions. The values of A are given in each panel. 


V. CONCLUSIONS 

We have analyzed the unitary time evolution of fi¬ 
delity and information entropy in an isolated interacting 
many-particle quantum (fermionic and bosonic) system 
followed by a random interaction quench as a function 
of many-particle spin. We obtain analytical results for 
time dynamics of fidelity and entropy in the framework 
of EE(l-|-2)-5' which are in good agreement with the nu¬ 
merical examples considered in the paper. We showed 
that the behavior of the spectral variances with spin ex¬ 
plains the spin dependence of fidelity and entropy. There 
is also a dependence on the quench strength A which is 
explored in all details in the paper with respect to the 
three chaos markers (Ac, Xe, Xd) generated by EE(l-|-2)-S'. 
It is important to note that BW to Gaussian transition in 
emission intensities (fidelity decay) has been studied ex¬ 
perimentally in [41] and the results are explained by the 
widths of the energy distributions. We have also derived 
approximate formulas for the saturation value S{oo) and 


saturation time tgat for the information entropy. In 
summary, an isolated interacting many-body fermionic 
(bosonic) quantum system exhibits delay in thermal- 
ization with increasing (decreasing) many-particle spin 
S{F). Our results may contribute to formulation of a 
general statistical description of non-equilibrium quan¬ 
tum systems and in designing quantum computers. 

Several approximations are used to derive the formu¬ 
las for fidelity and entropy as exact formulation is non¬ 
trivial. In future, one should attempt to solve EE(l-|-2)-5' 
[also EE(l-|-2)] from the first principles. It will also be 
interesting to compare the formulas for tsat and other 
quantities given in this paper with the results of realis¬ 
tic many-body calculations such as those presented re¬ 
cently by Lode et al. [55]. Another important aspect in 
relaxation of quasi-integral systems is prethermalization. 
Recently prethermalization has been experimentally con¬ 
firmed for a one dimensional degenerate Bose gas [56, 57]. 
Therefore, it will be interesting to explore prethermaliza¬ 
tion using EE(I-|-2)-S' formalism. 
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FIG. 6. Plot of S{t) as function of time t for different F spin values for BEGOE(l+2)-F’. The points, marked by red ’+’ 
{F = 2), blue ’open square’ {F = 4) and black ’open circle’ {F = 5) correspond to numerical results. Similarly, the continuous 
curves represent the corresponding theoretical predictions. The values of A are given in each panel. 


The generic features explored in this paper consti¬ 
tute an important step towards a complete description 
of unitary time evolution of quantum systems of inter¬ 
acting particles quenched far from equilibrium. This for¬ 
mulation gives an overall picture of relaxation of com¬ 
plex quantum systems in the absence of complete knowl¬ 
edge about system dynamics [which is generally the case 
for complex physical (such as nuclear, mesoscopic and 
atomic) systems]. Moreover, results of the present anal¬ 
ysis may be useful in the study of the stability of a quan¬ 
tum computer against quantum chaos [58] and to ad¬ 
dress Loschmidt echoes in many-particle quantum sys¬ 
tems [46]. 

EE(l-|-2)-S' analysis of time evolution of few-body ob¬ 
servables is in principle possible and it will be consid¬ 
ered in a future publication. We also plan to utilize BE- 
GOE(l-|-2) ensembles with spin one degrees of freedom to 
study non-equilibrium dynamics of isolated finite quan¬ 
tum systems. A method for constructing this ensemble 
is given in [59] although this ensemble is more challeng¬ 
ing, computationally as well as analytically. Attempts 


will also be made to obtain an analytical understanding 
of the significance and magnitude of the parameter n in¬ 
troduced in Section 3.2. 

After completing the present work, we became aware 
of a recent paper [60]. Here, EGOE(l) generated by a 
random one-body hamiltonian (with both diagonal and 
off-diagonal matrix elements) is used to show analytically 
that the corresponding fermionic systems exhibit ETH. 
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APPENDIX 


We give the formulas for the variance propagators m, S) and Q(n, m, S) for EG0E(l+2)-S' and BEGOE(l+2)- 
F derived in [37] and [38] respectively using trace propagation method: 


P(0, m, S) = 


+ 


1 


n{n + i)/2 

'n + 2 




0(0- l)/2 [0 + 1 


Qi({l2}:m,5) + 


O^ + O + 2 
02+0 




Qi({2}:m,5) = 
Q2({2}:m,5) = 

Q^({1^} :m,S) = 


16(0-2) 
+ 80(m — 1)(0 — 2 to + 4)5^] 
O 


(0 + l)iP°(m,5)/16] [m“(m + 2)/2 + 52] , 
0(0 + 3)P°(m, 5)/32] + 1) - 5^] , 

[{n + 2)V\m,S)V‘^{m,S) 


Qmi^}:m,S) = 


[(30^-70 + 6)(5^)" 


8(0 - 2) 

+3m(m — 2)rrF{nrF — 1)(0 + 1)(0 + 2)/4 
+ - 3)(0 + 2) + 0(0 - 1)(0 + 1)(0 + 6)}] , 

P°(to, S') = [TO(m + 2) — 45^] , P^(m, S) = [3TO(m — 2) + 4S2] , 
p2(m,S) = 3m'"(m-2)/2-S2 , m'" = (o - = S(S + 1) . 


(26) 


g(0, m, E) = 1)(0 - 2(-l)^)(0 + 2) n''=^’f {m, F) 

7?.=i,/=o. ^ [(m + 2)mV2-(E^)] X^jm^F) 

^ ^ 8(0-2)(0-1)0(0+ 1) 

80(m-l)(0 + 2m-4)(F2) + (0-2)A’2(m,F)A’i(m,F) 

8 ( 0 - 1 ) 0 ( 0 + 1)(0 + 2)2 

./= 0 (^^ E) = [m*(m* - 1) - (E^)] A’0(m, E)/[80(0 + 1)] , 

./=i(^^ E) = I [(E2)] ^ (302 7^ 6)/2 + 3m(m - 2)TO*(m* + l)x 

(O- 1)(0- 2)/8+ [(e 2) /2] [(50 + 3)(0- 2)TOm* + 0(0- 1)(0 + 1)(0- 6)]}/ 
[(O- l)0(0 + 2)(0 + 3)] ; 

A’°(m, E) = [m(TO + 2) — 4 (E^M , ■l:’^(m, E) = [3m(m — 2) + 4 (E^])] , 

F) = 3(m - 2)m*/2 + (e 2) , m* = O + m/2 , (E^) = E(E + 1) . 


(27) 
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